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ABSTRACT 

Here we report the results of an investigation into the effects of ionizing ra- 
diation on commercial-off-the-shelf InGaAs and Si photodiodes. The photodi- 
odes were exposed to 30, 52, and 98 MeV protons with fluences ranging from 
10^ — 5 X lO^-*^ protons/cm^ at the Indiana University Cyclotron Facility. We 
tested the photodiodes for changes to their dark current and their relative re- 
sponsivity as a function of wavelength. The Si photodiodes showed increasing 
damage to their responsivity with increasing fluence; the InGaAs photodiodes 
showed significantly increased dark current as the fiuence increased. In addition, 
we monitored the absolute responsivity of the InGaAs photodiodes over their 
entire bandpass. Our measurements showed no evidence for broadband degrada- 
tion or graying of the response at the fluences tested. All measurements in this 
investigation were made relative to detectors traceable to NIST standards. 

Subject headings: Astronomical Instrumentation 
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Introduction 



The accurate calibratio n of optica 



objective of SNAP/JDEM 



WFIRST fiBlandford et al. 



( Levi et al- 



and near infrared photome try was a primary 



2010 



Aldering et al. 



20021 ) . since renamed 



201l[ ). a space-based mission to study dark energy. In 



particular, systematics-limited measurements of the dark energy parameters based on 
SNIa observations required broadband filter measurements with a top-down error budget 



of 2% color error and a 1.5% in-band error (IMostek 



20071 ). To achieve these ambitious 



goals SNAP/JDEM considered a multi-technique approach to focal plane calibration that 
included accurate stellar photometry and an onboard flat fielding illumination system. For 
the low-frequency spatial flats (L-flats) that characterize large scale, many-pixel variations 
across the focal plane, the technique favored was to dither large ensembles of stars across 



the f o cal plane and to measure the stars' offsets from the ensemble mean (Ivan der Marel 



2003 



Mostek 



20071 ). On the other hand, high-frequency spatial fiats (S-fiats), or small 



scale few-pixel variations across the focal plane, are typically monitored by a flat-fielding 
illumination system. For SNAP/JDEM, one version of this onboard calibration light system 
included filament lamps and/or pulsed LEDs, light sources that were to be driven at 
constant current. Since the accuracy of such high-frequency fiat fielding systems depends 
strongly on the capacity to compensate for variations in the illuminators, stable and 
precisely calibrated photodiode detectors were included to monitor the light sources. 

Modern calibration strategies to achieve fiat fielding goals as challenging as those 
required for the SNAP/JDEM high frequency flats are typically based upon monitoring 
photodiodes calibrated to the state-of-the-art by NIST or equivalent national standards 
laboratories. Since SNAP / JDEM was planned as a multiyear mission at the L2 Lagrange 
point, approximately 1.5 x 10^ km from Earth, we studied these monitoring photodiodes 
for their ability to survive in the radiation environment expected. As described by 
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Dawson et al. I ( 120081 ). solar protons dominate the radiation environment at L2. For 
instance, in a six year mission with an experiment carrying a minimum Aluminum 
equivalent shielding of 9.0 mm with an average shielding thickness of 38 mm, the focal 
plane monitoring photodiodes would be exposed to fiuences of order T = 10^ — lO^'' cm~^ 
in the range 10 - 100 MeV. 

In this investigation, we describe radiation hardness studies of several commercial 
off-the-shelf (COTS) InGaAs photodiodes and one Si photodiode irradiated with protons 
at energies 30, 52, and 98 MeV, and with fiuences up to J-" = 5 x 10^^ cm~^ at the 
Indiana University Cyclotron Facility (lUCF). In our detector-based calibration scheme for 
SNAP/JDEM, we proposed InGaAs photodiodes to monitor the NIR fiat-field illuminators 
and Si photodiodes to monitor the optical illuminators. The intent was to have all the 
photodiodes calibrated by NIST. Here we examine both the relative spectral response 
and the overall broadband response of InGaAs and Si photodiodes as a function of the 
irradiating proton energy and fiuence. Alt hough there have bee n radiation hardness 



studies of Si photodiodes, including CCDs (Dawson et al. 1 12008 



), there is surprisingly 



little information on InGaAs photodiodes. As proposed space missions push into the NIR 
(EUCLID, WFIRST), InGaAs photodiodes are hkely to find wider use in space-based 
experiments. 



2. Photodiodes and Radiation Effects Research Program 



In this section we describe the photodiodes tested and the facility where the radiation 
exposures were made. 
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2.1. Photodiodes 

The photodiodes tested in these studies were obtained from five optoelectronic 
manufacturers. All were COTS devices and are listed in Table [TJ The InGaAs photodiodes 
all have a responsivity cutoff at 1700 nm and the Si photodiodes have the typical 
responsivity cutoff at 1000 nm. When possible, the photodiodes were obtained in TO-5 
packaging, which has been shown to be robust to the shake and heat cycle testing required 
to qualify for space ffightQ. Only the InGaAs photodiodes from Advanced Photonics were 
unavailable with this packaging as a GOTS device. Since our plan was to expose the InGaAs 
photodiodes to protons at three energies, we obtained three InGaAs photodiodes from each 
manufacturer. We also obtained an additional PerkinElmer InGaAs photodiode that was 
not exposed to radiation and which was used as a control to determine systematic errors. 
We used 7 OSI Si photodiodes in this study, one for each ffuence we planned on using. We 
only needed 7 Si photodiodes because they were irradiated at only one energy. 

The spectral response measurements of the test photodiodes were made with respect to 
stable reference photodiodes whose response was calibrated by NIST. The NIST-calibrated 
reference InGaAs and Si photodiodes are listed in Table [2J 



2.2. Radiation Effects Research Program 



The photodiodes were irradiated at the Radiation Effects Research Program (RERP) 



facility at lUCF 



The lUCF cyclotro n produces 200 MeV protons with dosimetry better 



than 10% (jvon Przewoski et al.ll2005l ). The energies chosen for the tests, 30, 52 and 98 



MeV, were selected to span a reasonable range of the particle environment expected at L2 



^ http : / / www. osioptoelectronics. no / custom-oem- solutions / Space . asp 
^http : / / www. iucf . Indiana, edu / rerp / 
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( iDawson et al. II2008I ) from the range of energies available at RERP. Energies lower than 30 
MeV are not available. RERP has both wide band and narrowband beam configurations. 
The wide band beam degrades the energy of the 200 MeV protons to the program energy 
using a thick copper plate upstream of the target photodiodes. In the narrowband beam 
configuration, the protons first pass through a magnetic spectrometer that selects out 
a narrow range of momenta and then these selected protons pass through a beryllium 
degrader to obtain the program energy. The advantage of the narrowband beam is that 
the radiation hardness studies can focus on specific damage mechanisms. The advantage of 
the wide band beam is that the proton flux is two orders of magnitude greater than the 
narrowband beam, thereby reducing exposure times. At 45 MeV the energy resolution in 
the narrow band beam is approximately 3.3 MeV (FWHM); for the broadband beam, the 
energy resolution is approximately 13.5 MeV (FWHM). 

For fluences J-" < 5 x 10^ cm~^, we used the narrowband beam and exposure times were 
a few minutes long. For fluences J-" > 5 x 10^ cm~^, we used the wide band beam to keep 
the exposure times reasonably short. Exposure information is given in Table [3l During 
radiation exposure, all of the target photodiode's pins were connected and grounded. 



3. Dark Current and Relative Spectral Response 

These investigations were primarily aimed at understanding whether the spectral 
response of photodiodes degrades in response to the ionizing radiation environment expected 
at L2. 
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3.1. Indiana Relative Responsivity Measurement Apparatus 

We developed the Indiana Relative Responsivity Measurement Apparatus (IRRMA) 
to measure dark current and relative spectral response as a function of wavelength; the 
apparatus is shown in Figure [TJ The apparatus consists of a radiance-controlled QTH 
lamp source powered by a const ant- current power supply, a monochromator that feeds into 
an integrating sphere, and a dual-channel picoammeter that measures the output current 
of the test photodiode and a NIST-calibrated reference photodiode simultaneously. The 
NIST-calibrated photodiode is temperature controlled. The light source (la) is a 100 W 
incandescent QTH lamp in a Newport Photomax housing. The lamp is powered by a 
radiometric power supply (lb) controlled by a Newport digital exposure controller. The 
Newport monochromator (2) has two gratings that span the wavelength range investigated 
here: 500-1000 nm for Si photodiodes and 1000-1600 nm for InGaAs photodiodes. The 
monochromator slits produce a 4 nm bandpass at the output. There is a shutter at the input 
of the monochromator that provides a measurement of the dark current when closed. The 
output of the 4" Labsphere integrating sphere (3) projects fiat illumination onto both the 
test photodiode and the NIST-calibrated reference photodiode. The two photodiodes are 
mounted in the dark box (4) where they simultaneously view the light from the integrating 
sphere. A Keithley dual-channel picoammeter (5) measures the two photodiodes with zero 
external bias. However, to make the current measurement, the picoammeter by design 
applies a measured 60 ± 5fiV reverse bias to the photodiode which in turn generates a few 
picoamps of dark current. 

The picoammeter sampled the output of the photodiodes at 6 Hz at every 2 nm 
wavelength step. To reduce read noise, 30 measurements were taken at each wavelength 
step and averaged. The dark current was measured every 10th wavelength step, or at 20 nm 
intervals. The dark current was subtracted from the signal at each 2 nm step by computing 
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a linear interpolation of the dark current over each 20 nm measurement interval. For the 
results described in §3.51 we made two full wavelength scans of each photodiode after an 
exposure at lUCF and we averaged the dark subtracted measurements. The dark current 
results we report below are the mean of the dark current measurements for both scans. 



3.2. Experiment Design 

After each radiation exposure, we measured the spectral response of the test photodiodes 
in IRRMA concurrently with respect to the stable NIST-calibrated reference photodiode. 
The reference photodiodes were calibrated by NIST with standard spectroradiometric 
detector calibration servicejf]. After each radiation exposure we made measurements of the 
test photodiode at each wavelength step and compared these to the reference photodiode 
at the same wavelength. To minimize systematic errors, we used the fractional change in 
the ratio of the responses of the test photodiode to the reference photodiode as the test 
statistic. 

The complication with reducing systematic errors in this way, however, is that 
differences between the response of the test and reference photodiodes can mask or 
exaggerate real spectral changes due to radiation exposure. If the spectral response of the 
reference photodiode at a particular wavelength is large compared with the test photodiode, 
for instance, significant spectral changes in the test diode would not change the ratio 
appreciably, thereby masking real spectral changes. On the other hand, if the spectral 
response of the reference diode is small compared with the test diode, small spectral 
changes in the test diode would lead to large changes in the ratio. To account for different 
spectral responses of the test and reference photodiodes, we normalized the ratios with 



^http: / / www.nist.gOv/calibrations/spectroradiometric.cfm#39075S 
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respect to those measured at Xnorm = 600 nm for Si photodiodes or at Xnorm = 1200 nm 
for InGaAs photodiodes. By normahzing the ratio in this way, we clearly only measured 
relative changes in the spectral response of the test photodiodes. 

To be quantitative, we defined AR\(t) as our test statistic. Let I\(t) = photocurrent 
response of the test photodiode to the QTH lamp at wavelength A measured in the 
A A = 4 nm monochromator bandpass at some time t during the experiment and let D{t) 
= the interpolated dark current at time t. Similarly, let A^a = photocurrent response of 
the NIST-calibrated reference photodiode to the QTH lamp at wavelength A and d = its 
interpolated dark current during the scan. Aside from variations in the lamp output, A^^ 
and d are assumed to have no additional time-dependent behavior. Then the normalized 
response of the test photodiode relative to the NIST calibrated photodiode at wavelength A 
at time t, Rx(t), is given by 



Rx{t) 



h{t)-D{t) 
N.-d 



hit) - D{t) 
N.-d 



'V. 



— ^norn 



where only the bracketed expression in the denominator is evaluated at the wavelength 
Xnorm- Wc quantified changes in photodiode response as a result of radiation exposure with 
the fractional change in Rx{t), 



ARxit) = Rx{t)/Rx{0)-l. 



(2) 



3.3. Systematic Uncertainties 

Several systematic effects in the apparatus that could introduce uncertainties that 
mask real changes in the spectral response of the test photodiodes cancel in the ratio 
defining Rx{t) in eq.([T]). For instance, QTH lamps have rated lifetimes of 50 hr and needed 
to be replaced on occasion during the course of these investigations. Since both the test 
photodiode and the NIST reference photodiode see the same lamp light, these systematics 
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cancel in Rx{t). One systematic effect not canceled in this way are variations in the 
test photodiode responsivity introduced by measurements in our laboratory, which is not 
equipped to maintain a constant temperature for the test photodiode. Since the responsivity 
of photodiodes as a function of temperature can differ depending on the manufacturer, 
this systematic effect would not be canceled in Rx{t). We tested for this systematic effect 
by measuring ARx for an unexposed control PerkinElmer InGaAs photodiode at 13. 4° C 
and 17.5°C, a temperature range greater than we encountered during our measurements. 
In these measurements, we attached a TEC to the back of the photodiode package and 
used an Omega controller to vary the temperature. The results are shown in Fig. |51 which 
shows clearly that negligible systematic variations of < 0.5% are introduced by temperature 
variations over this range. 

Before radiation exposure, we measured the dark current for all of the test photodiodes. 
The results of the measurements for the InGaAs photodiodes are given in Table |5l for the Si 
photodiodes, the dark currents are given in Table [61 Except for the Fermionics photodiodes, 
the baseline dark currents vary by up to a factor of ~50%; the Fermionics photodiodes show 
considerably greater variation. Again excluding Fermionics photodiodes, the dark currents 
vary by a factor of ~2 from manufacturer to manufacturer; the Fermionics photodiodes are 
by far the noisiest. 

We determined the systematic uncertainties in measurements of the dark current by 
repeatedly sampling the dark current of the control PerkinElmer InGaAs photodiode and 
assuming the dark current of this photodiode did not change. These measurements are 
shown in Figure |3l In the left panel, the figure shows typical variations in the dark current 
during a single full- wavelength scan that takes approximately one hour. The RMS/mean for 
these measurements is ~2%. The right panel shows measurements of the dark current over 
the duration of the experiment (~600 days) that only includes measurements of the control 
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PerkinElmer photodiode made at the same temperature. For these measurements, the 
RMS/mean in the dark current is ~12%. Adding the short term and long term variations 
in quadrature gives the systematic error in the dark current measurements of ~12%. 

3.4. Changes in Relative Spectral Response due to Radiation Exposure 

After each exposure, we used IRRMA to test our photodiodes for changes in relative 
spectral response. We considered three possibilities. First, it is possible that radiation 
exposure did not affect the relative spectral response. In that case, ARx{t) will be 
consistent with zero within measurement errors. However, if the relative spectral response 
did change, then eq.(Il]) shows that there are two different ways that radiation could 
affect the photodiode. In one, radiation exposure damages the spectral response I\ of the 
photodiode, resulting in a wavelength-dependent change in ARx- In the other, radiation 
exposure drives an increase in the dark current D, resulting in a wavelength- independent 
change in ARx. We have developed a simple statistic to differentiate among these three 
possibilities, 



where the sum is over the n discrete monochromator measurements in the wavelength 
range (Aq, A/), Tlx is the model for the radiation damage, and a is the RMS error in the 
measurement of ARx in our apparatus. 

Central to this analysis is the determination of a, since it sets the scale that 
differentiates between real physical effects and measurement error. We determined a for 
each photodiode individually from the initial ARx{0) measurements taken prior to radiation 
exposure. Since we have three instances of each photodiode, one for each proton energy, 
we constructed two ARx data sets - one in which the 52 MeV data were compared to the 




(3) 
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30 MeV data in eq.([2]) and one in which the 98 MeV data were compared to the 30 MeV 
data. Assuming the three instances of the the photodiodes are identical, we used the RMS 
of these two distributions to estimate of the a in eq.(|n]). Fig. H] shows the distributions of 
"identical" photodiodes for those manufactured by PerkinElmer. The value of a for the 
PerkinElmer photodiode, as well as those for the others in Table [H are given in Table HI 

We then used these a's to test against the three simple models for radiation damage 
with eq.(|3]). As usual, our strategy looks for /NDF ~ 1 for a successful model fit. For 
a photodiode that does not change, 1Z\ = 0. If x^/^df is significantly greater than 1 for 
the TZx = model, we assumed there was a change in the photodiode response and then 
tested the data against a model in which ARx(t) changes with A. For this model we set 7l\ 
equal to a fifth order polynomial fit to the ARx{t) data at each fiuence. If /NDF ~ 1 
for this model at each fiuence - that is, the ARx{t) data can be well described by a 
wavelength-dependent model within the measurement errors - we concluded that radiation 
exposure damages the spectral response Ix of the photodiode. If x^/NDF is significantly 
greater than 1 for this model, then we concluded that the changes in ARx{t) are due to 
wavelength-independent effects, like an increase dark current D as the fiuence increases. In 
this case, we looked for a correlation between increases in /ndf and dark current D{t). 

We demonstrate our analysis approach with an example. Fig. [5] shows ARx{t) for 
the OSI Si photodiodes irradiated at 30 MeV and fiuences between 5 x 10^ - 5 x 10^^ 
protons/cm^. There are clear and significant spectral response changes evident that increase 
with exposure. The left panel of Fig [H] shows x^/NDF for the TZx = model as a function of 
fiuence. As expected, this model fit is not a good one. The middle panel shows x^/NDF for 
the wavelength-dependent model of radiation damage. Since /NDF < 1 for this model, 
it is clear that radiation exposure damages the spectral response Ix of the Si photodiodes, 
with increasing damage as the fiuence increases. The Si photodiodes started off with less 
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than ~ 0.1 pA of dark current, as shown in Table [6l The right panel shows the variations in 
dark current with radiation exposure. Changes are well in excess of our ~ 12% measurement 
uncertainties. Since the photocurrents dominate the dark currents at all fluences, the dark 
currents do not affect our measurements of relative responsivity. 

The measured dark current after irradiation did not change in such a way to affect our 
ability to measure changes in the relative responsivity as a function of wavelength since our 
photocurrents are much greater than the measured dark currents after irradiation. 

3.5. Results for InGaAs Photodiodes 

We used the RERP facility to irradiate the InGaAs photodiodes in Table [T] according to 
the exposure history in Table [31 Before any radiation exposure, we scanned each photodiode 
and measured its dark current to determine -Ra(O) in eq.(l2]). The dark currents from these 
baseline measurements are given in Table O 

After each exposure we determined ARx(t) by rescanning the photodiodes and 
remeasuring the dark current. Fig. [7] shows ARx{t) for the PerkinElmer photodiode after 
exposure to 30 MeV protons at J-" = 10* and 5 x 10^^ protons/cm^. There appears to be 
little evidence for damage at 10* protons/cm^; there is, however, significantly more damage 
at 5 X 10^^ protons/cm^. The behavior seen at 30 MeV for this PerkinElmer photodiode is 
qualitatively similar to that seen in all the InGaAs photodiodes we investigated. 

We analyzed the ARx(t) data for the InGaAs photodiodes with eq.([3]) as was done 
for the OSI Si photodiodes. We first describe in detail the results from the PerkinElmer 
photodiodes. Since the behaviors of the remaining InGaAs photodiodes are similar, we give 
less detailed descriptions for them. Fig. [8] shows the results for the PerkinElmer photodiodes. 
The left panel of Fig [8] shows x^/NDF for the TZx = model as a function of fluence for the 
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three proton irradiation energies. ARx{t) shows httle evidence for damage below J-' = 10^° 
protons/cm^; above = 10^^ protons/cm^, however, the value of x^/NDF increases 
significantly with fluence, a result that implies increasing damage to the photodiode with 
radiation exposure. The middle panel shows x^/NDF for the wavelength-dependent TZx 
model, where AR\{t) is fit with a fifth order polynomial. The wavelength-dependent model 
clearly does not account for the increasing radiation damage above 10^° protons/cm^, which 
implies that the degradation in ARx(t) is due to wavelength- independent effects. The 
right panel shows the (negative) dark current as a function of fluence for the three proton 
energies. The increases in dark current are well above the 12% systematic error in the 
measurements. The correlation between the rise in x^/NDF of ~ 3 orders of magnitude 
above 10^° protons/cm^ and the comparable rise in wavelength-independent dark current 
is apparent. The x^/NDF behavior in the middle panel can be qualitatively understood 
as follows. I\ is relatively insensitive to radiation damage at the fiuences investigated 
here. At low fluence, Ix{t) dominates D{t) and ARx{t) shows little evidence of radiation 
damage. As D{t) increases with fluence, D grows to dominate Ix and dark current adds 
increasingly large fluctuations to each monochromator wavelength measurement. These 
fluctuations are well in excess of a. Since the sensitivity of COTS InGaAs photodiodes falls 
at longer wavelengths, while magnitude of the fluctuations in dark current are independent 
of wavelength, the relative importance of the fluctuations in ARx{t) increases at longer 
wavelengths. This behavior in apparent in Fig. [71 

The results for the remaining photodiodes analyzed with eq.(l3]) are shown in Figs. [9]- 
[121 The behavior seen in Fig. [HI is qualitatively similar to what is seen for these photodiodes: 
radiation damage leads to significant increases in wavelength-independent dark current that 
affect photodiode response. More specifically, the PerkinElmer and OSI photodiodes show 
the greatest increases in dark current. The API and Hamamatsu photodiodes show smaller 
increases. The Fermionics photodiodes show the smallest relative changes. However, the 



- 15 - 



baseline dark currents for the Fermionics photodiodes are much greater than the others. 
The increase in dark current as a function of radiatioii expo sure has also been reported for 



InGaAs Avalanche Photodiodes (IBecker and Johnson 



20041 ) 



4. Absolute Response of InGaAs Photodiodes through 10 nm Narrow Band 

Filters 

Since our spectral response measurements for InGaAs photodiodes only determine 
relative changes in response as a function of wavelength, these measurements would not 
detect radiation damage that affects the overall photodiode response in a wavelength 
independent way - a so-called "graying" of the response. To evaluate whether graying has 
occurred, we developed a second apparatus, the Absolute Responsivity Apparatus (ARA), 
that measures absolute changes in the responsivity of the InGaAs photodiodes with a 
filament lamp viewed through a set of narrowband interference filters. 



4.1. Absolute Responsivity Apparatus 

The ARA consists of a 50 W QTH lamp, an automated shutter, a filter wheel containing 
seven narrowband, 1.25" diameter Oriel interference filters, and hardware to mount the 
photodiode under test. The QTH lamp is positioned 0.25 m from the photodiode and is 
powered by a Newport constant-current radiometric power supply. We used the automated 
shutter to monitor the dark current before and after each filter measurement. The seven 
narrowband filters in the filter wheel are given in Table [71 The photodiodes were read out 
100 times per filter with a Keithley dual-channel picoammeter at a sample rate of 6 Hz to 
reduce read noise. The 100 measurements were averaged to determine the photo-current. 

The ARA yields photo-currents that are 10^ — 10^ times greater than the dark current. 
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This large signal dominates the baseline dark currents, as well as variations in the read-out 
electronics that can add uncertainty to the measurements, so that the ARA is testing the 
overall photodiode response. 



4.2. Experiment Design 

Since the filters are 1.25" in diameter, only one photodiode at a time can be measured 
in the ARA. We therefore mounted the test and reference photodiodes on a sliding stage 
so that they could be moved in and out of the apparatus. The measurements were carried 
out with the following cadence: the reference photodiode in all filters, the test photodiode 
in all filters, and then the reference photodiode in all filters once again. We averaged the 
two reference photodiode measurements in each filter that bracket the test photodiode 
measurements. Again we used the ratio of the response of the test photodiode to the 
average response of the reference photodiode to look for damage due to radiation exposure. 

In analogy with eq.([T]), we define the absolute spectral response of the test photodiode 
relative to the NIST calibrated photodiode in bandpass AA at time t. A^x{t), is given by 

/AA(t) 



AAx{t) 



(4) 



where lAx(t) are measurements of the response of the test photodiode, < N^x > is the 
average of the reference photodiode measurements flanking the test photodiode, and AA 
is the filter bandpasses given in Table [3 The fractional change in the absolute spectral 
response is then 

AAAxit) = AAxit)/AAxiO) - 1. (5) 

In our narrowband filter experiments we compared photodiode measurements after six 
different exposures - 1 x 10^ protons/cm^, 5 x 10^ protons/cm^, . . ., 5 x 10^^ protons/cm^ 
and we used the 1 x 10^ protons/cm^ exposure as the Aaa(O). The narrowband measurement 
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program was initiated after the 1 x 10^ protons/cm^ exposure. 

4.3. Systematic Uncertainties 

We determined the systematic uncertainties introduced by the filter-testing apparatus 
by once again measuring the unexposed control PerkinElmer InGaAs photodiode. We made 
five back-to-back measurements of this photodiode (without power cycling the apparatus) 
through the filters listed in Table [7] on seven separate days from October 2009 to May 
2011. We combined the control photodiode measurements on separate days to establish the 
response expected for a photodiode that did not change as a result of radiation exposure. 
We then compared this baseline behavior with the irradiated photodiodes as a test of the 
null hypothesis - is the irradiated photodiode consistent with the hypothesis that it did not 
gray? 

To model the response of a InGaAs photodiode that does not gray, we constructed 
a set of simulated experiments in each filter with the data from the unexposed control 
PerkinElmer photodiode. For these simulations, we first created a histogram of all possible 
AAaa values that pair a single measurement on one day with a single measurement on any 
other day (e.g., measurement #1 on day 1 with measurement #4 on day 3). We did not 
pair measurements on the same day because they were made without power-cycling the 
apparatus, unlike the measurements of the test photodiodes made after each exposure. 

For any given filter there are 525 AA^x entries in the histogram. We then constructed 
simulated experiments by drawing five values of AAaa from this histogram and computing 
their mean and standard deviation. This procedure simulates the statistics we computed 
through each filter for each photodiode after the radiation exposure program was completed. 
(Using the first exposure as the reference, there are five AAaa values from the six 
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exposures.) In all we constructed 10,000 experiments through each filter for a total of 
70,000 simulated experiments. 

In Fig. [131 we plot the results of the simulated experiments. This figure histograms 
the number of standard deviations that the mean is displaced from zero, < AA^x > /o", 
for all 70,000 simulated experiments. Superposed on this distribution is a Gaussian, the 
distribution expected for our simulated experiments if they differed from one another as a 
result of measurement error. The mean of this distribution falls close to zero, as expected 
in this case. The width is narrower than 1 because we are averaging five values that are all 
drawn from a normal distribution. 

Fig. [13] shows the behavior expected from a photodiode whose response has suffered no 
graying damage. 

4.4. Results for InGaAs Photodiodes 

We measured the absolute spectral response of the InGaAs photodiodes in Table [1] 
using the ARA. Since we initiated this program after the photodiodes had been exposed to 
a fluence of J-" = 10^ protons/cm^, we used the measurements at 10^ protons/cm^ as the 
baseline Aaa(O) in eq.Q. We then remeasured the test photodiodes after the remaining 
exposures in Table [3] The results are shown in Fig. [141 This figure is to be compared 
with Fig. [131 which shows the behavior expected from a photodiode whose response has 
suffered no graying. The shape and width of the data distribution in Fig. [TH is reasonably 
consistent with Fig. [131 The peak, however, is offset from zero. We infer this offset is the 
result of making an explicit choice for Aj\x{0) in AA/^x{t). In the simulated experiments, 
all days were treated equally in computing v4aa(0) in eq.(|l]), a procedure that mitigates the 
bias introduced by an explicit choice for A^x{0). We tested this hypothesis by choosing a 
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different day for v4aa(0) in evaluating eq.© and we found tliat tlie peak moved significantly, 
consistent with the hypothesis. 

We infer from the results shown in Fig. [H] that there was no graying of the InGaAs 
photodiode response at the radiation doses investigated here. This suggests there would be 
negligible damage affecting the absolute spectral response of InGaAs photodiodes with the 
radiation exposure expected at L2. 

5. Summary 

In this investigation, we report on radiation hardness studies of several COTS InGaAs 
photodiodes and one COTS Si photodiode exposed to ionizing protons with energies of 30, 
52, and 98 MeV, at fiuences up to J" = 5 x 10^^ cm"^ at the RERP at the lUCF. 

We scanned the relative spectral response and measured the dark current with 
the IRRMA apparatus we developed for this investigation. We found that both the Si 
and InGaAs photodiodes experience radiation damage as the fluence increases. The Si 
photodiodes showed wavelength-dependent radiation damage, particularly at wavelengths 
longer than A = 700 nm, primarily as a result of damage to their responsivity. The InGaAs 
photodiodes, however, showed evidence for wavelength-independent damage as a result of 
significant increases in their dark current as the fluence increased. 

We used the ARA apparatus we developed to measure absolute changes in the 
responsivity of the InGaAs photodiodes. This investigation was designed to determine 
whether there was an overall graying of the response of the InGaAs photodiodes after 
radiation exposure. By comparing the ARA measurements to simulated experiments 
constructed from baseline measurements of the InGaAs photodiodes before radiation 
exposure, we found that the test photodiodes showed little evidence for graying of their 
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response. 
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Table 1. Photodiodes tested. 



Manufacturer 


Semiconductor 


Active Area 


Package 


r art 




Type 


Diameter 


Tvpe 


Number 


Advanced Photonix Inc 


InGaAs 


1.5mm 


TO-39 


SD 060-11-41-211 


Fermionics 


InGaAs 


2mm 


TO-5 


FD2000W 


Hamamatsu 


InGaAs 


2mm 


TO-5 


G8370-82 


OSI Optoelectronics 


InGaAs 


3mm 


TO-5 


FCI-InGaAS-3000 


OSI Optoelectronics 


Si 


3mm 


TO-5 


OSD15-0 


PerkinElmer 


InGaAs 


2mm 


TO-5 


C30642GH 



Table 2. NIST calibrated reference photodiodes. 



Manufacturer Semiconductor Active Area Package Part 
Type Diameter Type Number 



Hamamatsu InGaAs 3mm TO-8 G5851-23 

Hamamatsu Si 3.6mm TO-5 S1336-44BK 
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Table 3. RERP exposure history at 30, 52, and 98 MeV. 



Date 


Cumulative Fluence 
[protons/cm^] 


Beam 
Configuration 


5/15/09 


1 X 10^ 


Narrow 


0/UD/Uy 


X iU 


Narrow 


7/31/09 


1 X lO'^ 


Narrow 


11/11/09 


5 X lO'^ 


Narrow 


2/22/10 


1 X 10^° 


Wide 


6/18/10 


5 X 10^° 


Wide 


8/09/10 


1 X 10^^ 


Wide 


9/24/10 


5 X 10^1 


Wide 



Table 4. Value for a used in calculation for each photodiode manufacturer determined 
from the basehne measurements of each of the photodiodes relative to the other. 



Photodiode 


a 


API 


0.009 


Fermionics 


0.02 


Hamamatsu 


0.01 


OSI 


0.01 


OSI Si 


0.002 


PerkinElmer 


0.008 
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Table 5. Dark current measurements for InGaAs photodiodes prior to irradiation; a 
different set of pliotodiodes was used for the 3 program irradiation energies. 



Photodiode 


30 MeV 


52 MeV 


98 MeV 


Manufacturer 


Dark Current [pA] 


Dark Current [pA] 


Dark Current [pA] 


API 


-0.98 


-1.45 


-1.53 


Fermionics 


-11.55 


-45.81 


-18.72 


Hamamatsu 


-3.02 


-3.20 


-2.75 


OSI 


-1.93 


-2.50 


-2.89 


PerkinElmer 


-1.73 


-2.83 


-1.92 



Table 6. Dark current measurements for Si pliotodiodes prior to irradiation; a different 
photodiode was exposed to 30 MeV protons at the fluences given in the first column. 



Fluence 


Dark Current 


[protons/cm^] 


[pA] 


5 X 10^ 


-0.073 


1 X 10*^ 


-0.065 


5 X 10^ 


-0.064 


1 X 10^° 


-0.11 


5 X 10^° 


-0.076 


1 X 10^^ 


-0.084 


5 X 10^^ 


-0.078 
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Table 7. Narrowband filters in the filter-testing apparatus. 



•^center 


FWHM 




(nmj 


(nmj 


% 


700 


10 


50 


800 


10 


45 


900 


10 


45 


1000 


10 


40 


1050 


10 


40 


1200 


10 


35 


1550 


10 


30 



Table 8. Systematic uncertainties in AAaa(^) for narrowband filter measurements. 



Acenter Simulated 
(nm) Mean 

700 -0.0003 ± 0.004 

800 0.002 ± 0.005 

900 -0.001 ± 0.007 

1000 -0.001 ± 0.004 

1050 0.0004 ± 0.01 

1200 0.009 ± 0.02 

1550 -0.004 ± 0.008 




Fig. 1. — The IRRMA apparatus for the measurement of dark current and relative spectral 
response. 
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Fig. 2. — Variations in the response of a Perkin- Elmer InGaAS photodiode introduced by a 
4.1°C temperature variation from 13.4°C to 17.5°C. The systematic variations are less than 
0.5%. 
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Fig. 3. — Variations in the dark current for the unexposed control PerkinElmer InGaAs 
photodiode. left panel: variations in the dark current for a typical single full-wavelength 
scan. The RMS/mean for this scan is ~2%. right panel: variations in the dark current seen 
over the ~600 days of the experiment. For measurements made at the same temperature, 
the RMS/mean is ~12%. 
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Fig. 4. — The ARx distributions for the PerkinElmer photodiodes before radiation exposure. 
Both the 52 MeV and 98 MeV photodiodes were compared to the 30 MeV photodiode. The 
RMS of these distributions was used to determine a in Table HI 
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Fig. 5. — ARx{t) for the OSI Si photodiode irradiated at 30 MeV, plotted at seven fluences 
between 5 x 10^ and 5 x 10^^ protons/cm^. 
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Fig. 6. — Model fits to the radiation exposure data for the OSI Si photodiodes. left: x^/NDF 
for the TZx = model as a function of fiuence for 30 MeV proton irradiation. Since x^/NDF 
^ 1 this model fits the data poorly, implying radiation damage, middle: x^/NDF for 
the wavelength-dependent model of radiation damage. Since x^/NDF < 1, we conclude that 
radiation exposure damages the spectral response of the Si photodiode. right: The (negative) 
dark current as a function of fiuence. Changes are well in excess of our ~ 12% measurement 
uncertainties. Since the photocurrents dominate the dark currents at all fiuences, the dark 
currents do not strongly affect our measurements of relative responsivity. 
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Fig. 7. — ARx{t) for the PerkinElmer photodiode after exposure to 30 MeV protons at 
1 X 10® and 5 x 10^^ protons/cm^. This behavior is representative of the behavior seen for 
all test InGaAs photodiodes at these fiuences. 



- 34 - 




10' 10" lO'' 10'° 10" 10'^ 10' 10" 10° 10'° 10" 10" 10' 10" 10° 10'° 10" io" 

Fluence [protons/cm^] Fluenoe [protons/om^] Fluenoe [protons/om^] 



Fig. 8. — Model fits to the radiation exposure data for tlie PerkinElmer InGaAs pliotodi- 
odes. left: x^/NDF for tlie Tlx = model as a function of fluence for the three proton 
irradiation energies. Below 10^° protons/cm^ Ai?A(^) shows little evidence for damage; above 
10^'^ protons/cm^, the value of x^/NDF increases significantly with fluence, which implies 
increasing damage to the photodiode with radiation exposure, middle: x^/NDF for the 
wavelength-dependent TZx model, where AR\{t) is fit with a fifth order polynomial. The 
wavelength-dependent model clearly does not account for the increasing radiation damage 
above 10^° protons/cm^, right: The (negative) dark current as a function of fluence for the 
three proton energies. The increases in dark current at fluences greater than 10^ protons/cm^ 
are well above the 12% systematic error in the measurements. The correlation between the 
rise in x^/NDF above 10^" protons/cm^ and the comparable rise in wavelength-independent 
dark current is apparent. This correlation would account for the increase in x^/NDF. 
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Fig. 9. — Model fits to the radiation exposure data for the Fermionics InGaAs photodiodes; 
see Fig. |8l left x^/NDF for the TZx = model as a function of fiuence for the three proton 
irradiation energies, middle: x^/NDF for the wavelength-dependent Tlx niodel. right The 
(negative) dark current as a function of fluence for the three proton energies. 




Fig. 10. — Model fits to the radiation exposure data for the API InGaAs photodiodes; see 
Fig. [HI left x^/NDF for the TZ\ = model as a function of fiuence for the three proton 
irradiation energies, middle: x^/NDF for the wavelength-dependent Tlx niodel. right The 
(negative) dark current as a function of fiuence for the three proton energies. 
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Fig. 11. — Model fits to the radiation exposure data for tlie Hamamatsu InGaAs pliotodiodes; 
see Fig. |8l left x^/NDF for the TZx = model as a function of fluence for the three proton 
irradiation energies, middle: x^/NDF for the wavelength-dependent Tlx niodel. right The 
(negative) dark current as a function of fluence for the three proton energies. 




Fig. 12. — Model fits to the radiation exposure data for the OSI InGaAs photodiodes; see 
Fig. [HI left x^/NDF for the TZx = model as a function of fluence for the three proton 
irradiation energies, middle: x^/NDF for the wavelength-dependent TZx model, right The 
(negative) dark current as a function of fluence for the three proton energies. 
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Fig. 13. — The behavior expected for a photodiode whose bandpass response does not gray. 
Plotted are the number of standard deviations that the mean of five measurements of AA^x 
is displaced from zero, < A^Iaa > /c, for 70,000 simulated experiments. Superposed is a 
Gaussian with mean -0.05 and standard deviation 0.6. 
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Fig. 14. — The behavior of InGaAs test photodiodes. Plotted are the number of standard 
deviations that the mean of five measurements of AA^a is displaced from zero, < AAaa > 
/cr, for the total data set. The shape and width of this distribution in Fig. [Tl]is consistent 
with the simulated experiments Fig. [131 The offset of the peak from zero likely results from 
making an explicit choice for Aa\{0) in AAA\(t), as described in the text. 



